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We present Quantum Monte Carlo calculations that probe the paramagnet-ferromagnet phase 
transition in a two-dimensional Stoner Hamiltonian. With a screened Coulomb interaction we ob- 
serve a first order ferromagnetic transition for short screening lengths, and a second order transition 
with a longer screening length, accompanied by a rising critical interaction strength. Finally, we 
discuss the consequences of our results for an ultracold atomic gas with finite ranged interactions. 
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I. INTRODUCTION 

Layered systems that display magnetic correlations 
have emerged as an important testbed of strongly cor- 
related physics. The Stoner Hamiltonian represents 
the simplest possible metallic system that undergoes 
a ferromagnetic transition. Since a mean-field analy- 
sis [1] predicts a second order transition, in the vicin- 
ity of the low temperature transition Hertz [5] pre- 
dicted that quantum fluctuations would drive criti- 
cal behavior. However, it has recently been pre- 
dicted that quantum fluctuations are even stronger 
than envisaged by Hertz and drive phase reconstruc- 
tion through a first order ferromagnetic transition [TJ [3J- 
8 , consistent with the phase transition observed in the 
quasi two-dimensional systems Sr3_ x Ca a: Ru207 [9j HO]. 
Sr 2 Ru0 4 [II], La^Sr^RuCU [S], Ca 2 Ru0 4 [T3], and 
UGe2 [T4TU6| . The importance of the quantum fluctu- 
ations in driving this phase reconstruction motivates a 
careful theoretical analysis of the Stoner Hamiltonian. 
However, analytical studies [TJ[3HE] of the magnetic tran- 
sition rely on a perturbation theory in the interaction 
strength, whereas in reality the interactions are strong 
and the quantum fluctuations dominant. To probe the 
phase transition in the non-perturbative regime we per- 
form and present the first Quantum Monte Carlo (QMC) 
calculations of the two-dimensional itinerant ferromagnet 
with short ranged interactions. We complement this with 
the first analytical study of the ferromagnetic transition 
in two dimensions with a screened Coulomb interaction 
that allows us to be the first to calculate the dependence 
of the tricritical point temperature on screening length. 
The QMC calculations employ a fixed node approxima- 
tion, tempered by backflow corrections, so should com- 
plement and improve upon the accuracy of the analytical 
findings. 

An ultracold atomic gas could be an attractive alter- 
native realization of the Stoner Hamiltonian. Recent ex- 
periments on the three-dimensional system |17j . and po- 
laron systems [HI HH] have delivered some evidence for 
ferromagnetic ordering [2"0H2"2"] , though a competing loss 
process provides an alternative explanation 23, 24]. Ac- 
curate QMC calculations that pin the transition down 
should help guide future experiments that could realize a 
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FIG. 1. (Color online) The phase diagram for the system 
from the paramagnetic to ferromagnetic phase (shaded gray) . 
The red solid line shows the first order ferromagnetic tran- 
sition and the dashed line the second order transition with 
the tricritical point highlighted by the blue dot. The three 
lines denote boundaries with a screening length of kpb 2 = 0, 



Kb- 1 = 0.2, and kitf 



0.4. 



Hamiltonian analogous to the Stoner system [Tl [71 |2"5T[2T3] . 
The theoretical and experimental study of the idealized 
two-dimensional ultracold atomic gas also presents an op- 
portunity to shed light on high temperature supercon- 
ductivity where antifcrromagnetism competes with the 
d-wave superconducting phase pQ. 

In this paper we present QMC calculations to ana- 
lyze the paramagnet-ferromagnet transition in two di- 
mensions. To connect to the solid state the QMC cal- 
culations are performed first with the screened Coulomb 
inter-particle potential and later the results are compared 
to the square well potential. Both inter-particle poten- 
tials are characterized by a range parameter that we vary 
to gauge the consequences of screening in the solid state, 
and the interaction effective range in the cold atom gas. 
Combining our QMC results with a complementary an- 
alytical order-by-disorder approach allows us to derive 
the phase diagrams shown in Fig. [I] at finite tempera- 
tures. The phase diagram shows that the paramagnet- 
ferromagnet transition reverts from first to second order 
on increasing the screening length, and the correspond- 
ing critical interaction strength increases. The introduc- 
tion of a finite interaction range with kpb > increases 
the critical interaction strength, and lowers the tricritical 
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point temperature to be in line with experimental values. 
Finally, we adapt our formalism to assess the opportu- 
nity to observe ferromagnetic correlations in an ultracold 
atomic gas. 

II. FORMALISM 

To model the ferromagnet we focus on the modified 
Stoner Hamiltonian 

H = Y^ £pnp<T + JJ dr t dr ;ff (r t - r^h^n^ , (1) 

p,rx 

where e p = p 2 /2 is the dispersion, n prT is the fermion 
occupancy in momentum space, n ra in real space, spin 
o £ {t>4-}> an d we adopt atomic units with H = m = 1 
throughout. To simulate the inter-particle repulsion we 
use a screened Coulomb interaction parameterized as 
g{r) — ge~ r / h /2nbr acting between opposite spin par- 
ticles separated by a distance r. It is characterized by 
a screening radius b and interaction strength g To val- 
idate our results we also study the square potential of 
radius R, g(r) = gO(R — r)/irR 2 with Heaviside function 
9, whose interaction strength and screening parameter 
are related to the screened Coulomb interaction through 
a momentum-space expansion [30] . The interparticle po- 
tentials are defined so that they have the simple momen- 
tum space forms g q = g/\/l + b 2 q 2 and g(l — R 2 q 2 /8) 
respectively. This allows them to be used within our an- 
alytical formalism. These definitions mean that the fer- 
romagnetic transition in the mean-field approximation 
emerges at the critical interaction strength <?mf = 2tt. 
To evaluate the energy of the electrons we employ two 
complementary techniques, a Quantum Monte Carlo cal- 
culation limited only by a fixed node approximation, and 
an analytical evaluation of the free energy derived by 
a functional integral formalism. The mechanics of both 
methods are outlined below before we study the resulting 
phase diagrams. 

A. Quantum Monte Carlo formalism 

To seek and calculate the ground state of the Hamil- 
tonian we perform QMC simulations with the code 
CASINO [31] ■ This method optimizes a trial wave function 
at zero temperature, and finds the exact ground state 
subject to the nodal surface of the wave function be- 
ing fixed. The approach is a refinement of that used in 
previous studies of itinerant ferromagnetism [BJ I32H37| . 
We use a variational wavefunction ip = e~ J D that is 
a product of a Slater determinant, D, which takes full 
account of the fermion statistics and includes further 
electron-electron correlations through a Jastrow factor 
J. The QMC simulations comprise of two stages: firstly 
in variational Monte Carlo (VMC) the ground state en- 
ergy was minimized by varying the parameters in the Jas- 
trow factor; secondly diffusion Monte Carlo (DMC) starts 
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FIG. 2. (Color online) The error in energy AE for the k = 
(blue) and twist averaged (red) calculations with system size 
(N) for the screened Coulomb potential with b — 0, ff/flMF = 
0.5. The red dotted line shows the ~ 1/N 2 scaling. The green 
dashed line shows the energy resolution required to determine 
the order of the ferromagnetic transition. 



from the VMC wave function, and treats the Schrodinger 
equation as a diffusion equation to project out the exact 
ground state subject to a fixed node approximation. 

The Slater determinant D = det({V>kefc Ft , ^kefeFj.}) 
consists of plane-wave orbitals ^k(r) = exp(ik • r) whose 
momenta k satisfy periodic boundary conditions in the 
square simulation cell, and lie within the up/down spin 
Fermi surfaces /cf<t ■ With a square simulation cell, to en- 
sure that within VMC we have a real-valued wave func- 
tion and that the Fermi surface is circular, the num- 
ber of states must correspond to closed shells contain- 
ing iV CT = {1, 5, 9, 13, 21, 25, 29, 37, 45, 49} electrons, thus 
constraining us to discrete values of the magnetization. 
For computational efficiency we factorize the Slater de- 
terminant into up and down-spin components |38| . so 
D = det({i/'kGfc Ft }) det({-0kefc F4 })- Provided that the or- 
bitals of the minority spin state are the lowest energy 
orbitals of those in the majority spin state, [39] . this is 
the state with total spin s — s z = (JV-j- — N\)/2, 

The Jastrow factor, J({ri}), accounts for electron- 
electron correlations. It has the general form 

JV-l N 

J{{n})= E M M+P(ru), (2) 

where the summation over indices {i,j} cover all N = 
Nf+N± electrons, the electron separation is = r — rj, 
and Tij = |ry|. The Jastrow factor includes the polyno- 
mial expansion in electron-electron separation proposed 
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in Ref. 
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chosen so that it satisfies the Kato cusp conditions at 
rij = 0, is zero beyond the cutoff length L imposed by the 
Heaviside function 0, = 1/4 for equal spin electrons 
and Tij = 1/2 for opposite spins, and contains iV u = 8 
variational parameters {ai}. The cutoff length L was 
chosen to be the largest circle that could be inscribed in 
the simulation cell. The Jastrow factor also includes a 
plane- wave expansion 
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where the {G^} are the reciprocal lattice vectors of the 
simulation cell belonging to the Ath shell of vectors un- 
der the full symmetry group of the Bravais lattice, and 
the superscript "+" means that if is included in 

the sum then —Ga is excluded. The summation over 
A covers N p = 8 shells, with corresponding variational 
parameters {a^}. To broaden the freedom of the vari- 
ational wave function we also include backflow correc- 
tions [H]. These substitute the electron coordinates r, 
in the Slater determinant a new set of collective coordi- 
nates Xi = r,; where the backflow displacement 
of electron i is given by 
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where L is the cutoff length, N n — 7 the expansion order, 
and {cfc} the variational parameters. The inclusion of 
the backflow corrections allows the nodal surface of the 
wave function to shift and therefore relax the fixed node 
approximation. To seek the ground state the variational 
parameters of the trial wave function were numerically 
optimized within VMC by minimizing the total VMC 
energy 42 . The optimized VMC wave function was used 
as the trial state for the DMC calculation. 

The DMC method 38i simulates a population of walk- 
ers whose evolution is driven by the imaginary time 
Schrodinger equation to project out the ground state 
component of the VMC wave function. The walk is taken 
in a series of discrete time steps, with walkers branching 
or annihilating according to the local energy. The choice 
of time step, the control of the walker population, and 
the system size can each introduce errors into the final 
prediction of the ground state energy. We therefore now 
address how to minimize each source of uncertainty in 
turn. 

Time step: One can propagate forwards in time ex- 
actly by using Green's function Monte Carlo 



but unfortunately the method is computationally expen- 
sive |31[ 15 lj . Therefore here we employ an approximate 
Green's function that would become exact in the limit 
of short time-steps t for the walker evolution. How- 
ever, the computational effort required to achieve a given 
uncertainty in the prediction for the ground state en- 
ergy increases as 1/r, ruling out the use of infinitesimal 
time steps in practice. Therefore, where high accuracy 
is required, we use two different finite time steps {r^} 
and extrapolate to r = to obtain the ground-state en- 
ergy. Tests revealed that the time step error had entered 
the linear regime at r < 0.01 (at a Wigner-Seitz radius 
r s = 1). Here to minimize the uncertainty in the extrap- 
olate we followed the prescription of Lee [52] and ran 
simulations with time steps r = 0.01 and r/4 = 0.0025 
for relative durations of 1 : 8 respectively, and then fi- 
nally extrapolated to r = 0. 

Population control: To ensure that the population con- 
trol bias is negligible in all runs the target population 
exceeded 2000 |-jTJ[ . The number of equilibration steps 
discarded during the equilibration phase of each DMC 
calculation was set so that the root-mean-square distance 
diffused by a typical electron exceeded the simulation cell 
size. 

System size: The VMC and DMC simulations must 
be performed in a finite sized simulation cell that is peri- 
odically repeated to create an effectively infinite system. 
However, using the finite sized simulation cell introduces 
errors into the final prediction for the energy. The error 
in the estimated energy due to the finite sized simula- 
tion cell can be divided into three components: single 
particle kinetic energy, Hartree energy, and exchange- 
correlation energy. The error in the kinetic energy arises 
because of the approximation of the circular Fermi sur- 
face by the discrete set of k-vectors of the closed shells 
within the simulation cell. As the system is enlarged 
the k-vector grid becomes more fine resulting in abrupt 
changes in the kinetic energy. The use of a non-zero 
simulation cell Bloch wave vector k s that causes some 
k-vectors to lie outside of the Fermi surface, and others 
within can lead to a dramatic reduction in the finite size 
error, with the optimal being the Baldereschi point [54] . 
However, it is even better to take an average expecta- 
tion over all Bloch vectors k s within the first Brillouin 
zone [53]. We adopted the most optimal strategy which 
is to average over a uniform grid of twists (we use 1000 
points) centered on the Baldereschi point of the simula- 
tion cell Brillouin zone 53 . The Hartree energy is neg- 
ligible due to the short-ranged nature of the interaction 
whose maximum exponential decay length is 1/20 of 
the typical simulation cell size. The final contribution to 
the finite sized error is due to long-range two-body cor- 
rections. The change in kinetic energy can be encapsu- 
lated by the Chiesa-Holzmann-Martin-Ceperley approx- 
imation [55] , It has been shown that in two dimensions 
the scaling with long ranged interactions is ~ 7V~ 5 / 4 
whereas the with short ranged interaction such as our 
screened Coulomb interaction the scaling is expected to 
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be ~ A'' -2 [371 [S3]. We also note that it has been found 
that other expectation values such as magnetic suscepti- 
bility can suffer more from finite size effects, but here we 
focus only on the ground state energy [56] . 

To gauge the scale of the finite size error we ran tests 
on a paramagnetic system with the screened Coulomb 
interaction with 6 = and g/gwiF — 0.5 and vary the sys- 
tem size by changing the number of electrons N. Fig. [2] 
shows that twist averaging delivers an error in the energy 
(determined against a 1 /N 2 extrapolation to the infinite 
sized system at 1/N = 0) that is over an order of magni- 
tude smaller than that from taking just the single result 
calculated at k = 0. The scaling 1/N 2 appears sim- 
ilar to that found in a previous two-dimensional study 
with short-ranged interactions. [37] For the typical sys- 
tem size used 2/N < 0.02 the finite sized error is almost 
two orders of magnitudes smaller than the smallest en- 
ergy scale AE ss 0.012 that we will need to resolve the 
features of the phase transition (see Fig. [3ja,b)). For 
large magnetization the energy associated with the mi- 
nority spin species will have a finite size correction that 
scales as V-^minority but since that species now makes 
only a small contribution to the overall energy it is be- 
yond the order needed. Therefore twist averaging ensures 
that finite sized errors are inconsequential when analyz- 
ing the phase diagram. Finally we note that all of our 
twist-averaged energies are always slight over estimates 
of the true energy since the single-particle kinetic energy 
k 2 /2 is a convex function and the occupied k-space is a 
convex polyhedron [53j . 

Changing polarization: Previous studies of 3 He have 
established that Slater- Jastrow wave functions overesti- 
mate the unpolarized state energy, thus favoring the fer- 
romagnetic state. Including many-body backflow correc- 
tions can help reduce the energy of the unpolarized state 
to better align with experiment [57, 58]. Here we include 
two-atom backflow corrections, that for short ranged in- 
teractions have previously been found to have a rela- 
tively small impact in the ground state energy versus that 
for long-ranged interactions [37 , 59 , indicating that the 
QMC bias towards unpolarized states is less for shorter 
ranged interactions. However, we also note that an alter- 
native calculation on a lattice with a smaller fixed node 
bias indicates that the ferromagnetic transition is infinite 
order rather than first order. |60j 



turn partition function from a coherent state field integral 



B. Analytical formalism 



The itinerant ferromagnet has previously been ana- 
lyzed using a functional integral formalism. This will 
provide a useful complementary tool to study the ferro- 
magnetic transition. The formalism calculates the quan- 
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where the field if) describes a two component Fermi gas. 
We now decouple the quartic interaction term with a 
Hubbard-Stratonovich transformation into the full vec- 
tor magnetization <f> and the density channel p [7J. This 
is the simplest decoupling scheme that maintains rota- 
tional spin invariance and yields the correct Hartree- 
Fock equations [611 E2] • With the action now quadratic 
in the Fermionic degrees of freedom we integrate them 
out to recover the quantum partition function Z — 
J T><pT>pexp(— S) with the action given by 



S = Tr [<f>g(j)— pgp] — Trln (d T +p 2 /2 — p+gp)I— g<p-cr 



Here we have employed the operator form of the inter- 
particle potential g. We adopt M and po as the puta- 
tive saddle point values of the fields <fi and p, and then 
expand to quadratic order in the fluctuations in those 
fields about the assumed saddle point. This allows us to 
integrate out the fluctuations in the magnetization and 
density channels whilst keeping terms up to second order 
in the interaction strength [6l El [63] . This yields the free 
energy [7] 
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where the first term corresponds to the kinetic energy 
(with a Fermi distribution n p<J = 1/(1 + e^ p "/ T ), £ p(T = 
p 2 /2 — jit — agM, and T is the temperature), the second is 
the mean-field contribution of the interactions, and the 
third higher order interaction effects. The first two terms 
would be delivered by the standard Stoner mean-field 
theory, and the final term is attributed to fluctuation cor- 
rections [7J. We have opted, without loss of generality, to 
set the quantization axis along the direction of the mag- 
netization. This expression for the free energy remains a 
function of the magnetization and density, that can now 
be determined by minimizing the free energy with respect 
to M and p, thereby fulfilling our premise that these are 
the saddle point values. Since the Fermi distributions 
have a temperature dependence we can use the formal- 
ism to not only study T = 0, but also the phase behavior 
at finite temperature. The formalism applies not only 
at zero temperature, but also at finite temperature, thus 
allowing us to map out the entire phase diagram. 

The screened Coulomb interaction in momentum space 
is g(p) — gj y/l + b 2 p 2 . We have included the effect 
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(a) First order, fe|6 2 = (b) Second order, k^b 2 = 0.2 (a) Zero temperature phase behavior 
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FIG. 3. (Color online) (a,b) Energy bands as a function of 
magnetization M at the interaction strengths shown in (c,d) 
for the screened Coulomb potential (the lowest band is the 
weakest interaction strength). The red points highlight the 
minimum point in each energy band. The error bars are ap- 
proximately the line width. (c,d) Magnetization as a function 
of interaction strength across the transition. The left hand 
plots are at kpb 2 = 0, and the right-hand at kpb 2 = 0.2. 



of finite ranged interactions into the free energy fol- 
lowing the prescription in Ref. 63 , where it was ap- 
plied to a three-dimensional system. The mean-field 
term interaction strength is independent of momentum 
exchange so is left unaffected by the screening length 
b. In the fluctuation correction term the denominator 
means that dominant contributions to the momentum 
summation arise at |pi — P3I = \/2kF (6j [63] so that 
we can simply adopt this fixed value within the interac- 
tion g(pi — P3) 1 — ^ g(^kp), and therefore in the pres- 
ence of screening the fluctuation correction term is sim- 
ply rescaled by a factor of 1/(1 + 2k F b 2 ) 2 [63j . Here 



we use the definition kp — y'3TT 2 (n^ + nj.) where n a is 
the density of the electrons with spin a. Following this 
rescaling, the ground state magnetization can again be 
extracted by minimizing the free energy with respect to 
magnetization. 



III. RESULTS 

With both the Quantum Monte Carlo and analytical 
formalism in place we are well positioned to study the 
emergent phase diagram. We first focus on the screened 
Coulomb potential relevant for the solid state, before 
looking at the cold atom gas. 
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(b) Zero temperature first order transition 
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(c) Finite temperature phase behavior 
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FIG. 4. (Color online) (a) Phase boundary for the emer- 
gence of ferromagnetism with varying interaction strength 
(g/gMF, primary x-axis) and screening length (kpb 2 , primary 
y-axis) at zero temperature. The secondary x-axis and y- 
axis show the parameters relevant for cold atom gases. The 
solid line corresponds to a first order transition, and the 
dashed line a second order transition. The red line shows 
the boundary for the screened Coulomb interaction, and the 
blue line for the square potential. The magenta solid line 
shows the analytical result for the first order transition with 
a screened Coulomb interaction, (b) The magnetization at 
the first order transition as a function of Coulomb screening 
length calculated with QMC (red) and analytically (blue). 
The green line shows the minimum magnetization of a first 
order transition that can be resolved within QMC. (c) The 
shift in the phase boundary interaction strength, g — <?t=o 
from the screened Coulomb interaction zero temperature an- 
alytical phase boundary shown in ((a), magenta solid line) at 
finite temperatures T G {0, 0.1, 0.2, 0.3}Tf calculated using 
the analytical formalism. The red solid lines denote a first 
order transition and the dotted red line a second order transi- 
tion, and the green dots the tricritical points. The inset shows 
the tricritical point temperature with screening length. 
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A. Solid state 

We first study the screened Coulomb interaction rel- 
evant for the solid state. With a short ranged interac- 
tion k^b <C 1, the Fourier transform is momentum inde- 
pendent so we recover the contact interaction limit. We 
performed DMC calculations to determine the ground 
state energy for several different values of magnetiza- 
tion and interaction strength, plotting the energy bands 
in Fig. [3ja). The minimum in each band reveals the 
magnetization at that interaction strength. With rising 
interaction strength the magnetization jumps from zero 
to M fa 0.45, characteristic of a first order transition. 
This gives the magnetization with interaction strength 
curve in Fig.[3]jc) that demonstrates not only the first or- 
der transition, but how the magnetization subsequently 
grows into the fully polarized state. The prediction of 
a magnetic transition at g/gMF = 0.57 is in good agree- 
ment with the Eqn. Q prediction of a first order transi- 
tion at g/gMF = 0.51, that was also found in a previous 
analytical study [T]. 

We next repeat the procedure with a larger screening 
length klb 2 = 0.2. Fig. [3jb,d) shows that the magneti- 
zation grows smoothly with interaction strength, demon- 
strating a second order transition within the magnetiza- 
tion resolution of our QMC calculations. We summarize 
our results for several screening lengths in the phase dia- 
gram in Fig. |4^a). This demonstrates how upon increas- 
ing the screening range the critical interaction strength 
increases and the transition reverts from first to second 
order. This trend is in good agreement with the ana- 
lytical predictions from Eqn. (j7|, though the analytical 
formalism predicts that the transition remains first order 
due to the presence of a logarithmic divergence in the 
free energy F ~ \M\ 3 log(T) [HUE]. 

In order to probe the question of the order of the tran- 
sition more closely, in Fig. [4jb) we show the magneti- 
zation formed at the first order transition as a function 
of the screening length kpb 2 , predicted by both QMC 
and from the analytical formalism. The analytical for- 
malism predicts a magnetization that falls rapidly with 
screening length because the screening reducing the fluc- 
tuation correction term that was responsible for driving 
the first order transition at a lower interaction strength, 
but the transition remains first order due to the loga- 
rithmic divergence in the free energy. The magnetization 
following a first order transition with the contact inter- 
action is predicted to be smaller by QMC, and drops to 
zero at k\b 2 = 0.2, at which point the transition be- 
comes second order. However, since the QMC method 
is able to sample the magnetization only at discrete val- 
ues, M € {0,0.07,0.19,0.28,...}, the minimum magneti- 
zation that can be formed in a first order transition is 
M min — 0.19. Due to the rapid decay in the magne- 
tization formed following the first order transition with 
screening length, we cannot resolve whether QMC pre- 
dicts a first or second order transition at k\b 2 > 0.1, how- 



predicted at k^b 2 fa 0.1 by QMC and k^b 2 w 0.18 from 
analytics is consistent. 

The DMC results have good qualitative agreement 
with our analytical formalism at T = 0. This gives us 
confidence to use the analytics to study the phase be- 
havior at finite temperature. Fig. [I] and Fig. |4|c) show 
that the interaction strength of the transition increases 
with temperature as the fluctuation correction term is de- 
pressed pQ. The trend remains the same with increasing 
screening length. However, with increasing temperature 
the transition reverts from first to second order, and the 
inset of Fig.Qc) reveals that the tricritical point temper- 
ature reduces markedly with increased screening length. 



B. Cold atom gas 

We now study the possible emergence of ferromagnetic 
correlations in an ultracold atomic gas. The electrons are 
simulated by a two-component gas of fermionic atoms, 
that are mapped onto pseudo up and down-spin species. 
The interactions between the atoms are controlled using 
a Feshbach resonance. This gives experimentalists ac- 
cess to not only a range of different interaction strengths 
but also the effective range of the interparticle poten- 
tial |18j . that can be either positive or negative depend- 
ing on both the elements and the Feshbach resonance 
used. In a cold atom gas the parameter used to char- 
acterize the interaction strength is the two-dimensional 
scattering length a2D that can be related to the inter- 
action strength used within the Coulomb interaction by 
ln(fcpa2D) = 5mf/27t<7, and the effective range is linked to 
the Coulomb interaction screening length by r c = —b 2 /a. 
However, this gives us access to only negative effective 
ranges r c < 0. To study a positive effective range we 
turn to the square potential U6(R — r) to model the 
interactions. Not only will this give allow us to study 
positive effective range, but it will also serve as a use- 
ful point of comparison to the results for the screened 
Coulomb interaction. The square potential has a two- 
dimensional scattering length a2D given by ln(fcpa2D) = 
(ffMF/87r 2 i?)/[l - tanh(x)/x] w hh X = RVU, and effec- 
tive range 
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ever the predicted crossing of the boundary M n 



0.19 



We now use DMC calculations to determine the ground 
state in the presence of the square well potential. This 
allows us to augment the phase diagram in Fig. [4] The 
critical interaction strength and the emergence of the sec- 
ond order transition are similar for both the square well 
potential and our previous study of the screened Coulomb 
interaction, though the square potential can be extended 
to effectively negative kpb 2 . This verifies the form of 
our phase diagram, and confirms that irrespective of the 
inter-particle potential our main result of the transition 
changing from first to second order with rising screening 
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length within the resolution of our simulations, accom- 
panied by a rise in the critical interaction strength is 
robust. 

The persistence of the ferromagnetic transition from 
negative into positive effective ranges kpr c makes the 
two-dimensional gas a tantalizing target in the search 
for ferromagnetism in a cold atom gas. Recent experi- 
ments have shown that atomic gases with a large nega- 
tive effective range [JH] have greatly suppressed two and 
three-body losses. Our predictions of the changing mag- 
netization in Fig. [3] and Fig. ^b) could be measured 
through phase contrast imaging that would reveal the 
growing polarization in the emerging domains |22| . Our 
accurate determination of the phase boundary through 
DMC calculations, complemented by analytical predic- 
tions, should guide future experiments. 

IV. DISCUSSION 

In this paper we have outlined Quantum Monte Carlo 
calculations on the two-dimensional itinerant ferromag- 
net. The calculations show that the system undergoes a 
paramagnet to ferromagnet phase transition. The tran- 
sition is first order for interactions with a short effective 
range, and within our magnetization resolution second 
order for interactions characterized by longer effective 



ranges. We have reconciled this with the damping ef- 
fect that long range interactions have on the quantum 
fluctuations. We have compared the DMC results to an 
analytical formalism and found good agreement, allow- 
ing us to then use the analytic formalism to extend the 
phase diagram to finite temperatures. Higher tempera- 
tures suppress the first order behavior leading to a tri- 
critical point at T = 0.3Tp, whose temperature reduces 
with increasing screening length. 

With typical solid state experimental systems having 
a screening length 0.5 < kpb < 1 [TTJIM], the first order 
transition should still be visible but with the fluctuation 
correction term tempered, and thus the associated tri- 
critical temperature will be reduced. This could reduce 
the tricritical point seen in theory 0.3Tp to the 0.02Tp 
seen in experiments on both two (and three) dimensional 
materials [9 - 16J and increase the visibility of exotic low 
temperature phases that are promoted by quantum criti- 
cality such as the spin spiral and p-wave superconductor. 
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